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The Naval Postgraduate School has conducted five research 
cruises in California coastal waters for which sufficient 
data was obtained to allow boundary layer mixing rates to be 
determined. These data have been processed to determine the 
mixing rates. The rates have been correlated with meteoro- 
logical conditions and geographical location and average 
values for use in air pollution models have been determined. 
A simplified method for calculating the mixing rate from mean 
meteorological parameters is presented. 
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I. Introduction 

The rate at which air is mixed throughout the atmospheric 
boundary layer is one of the determinating factors in air quality. 
The mixing rate governs how rapidly pollutants from surface 
sources are transported upward, and thus is one of the factors 
that determines peak concentration in a given volume. Air 
pollution models divide a particular air shed area into numerous 
sections (cells) then solve for temporal variations of pollutant 
density within each cell using source strength and diffusion 
within the cell and flux across the cell boundaries as input 
parameters. Matching conditions at the boundaries leads to an 
overall solution to the problem. The transport across cell 
boundaries depends mainly on mean wind speed. The transport out 
of the top of a cell and vertical dispersion within a cell are due 
to turbulent mixing. 

Many of the model cells are over the water for air sheds along 
the California coast. In the past, there has been little data 
available that would allow over ocean mixing rates to be 
determined so estimates of the rate have been based on overland 
values. The purpose of the work reported here is to determine the 
needed mixing rates from available Naval Postgraduate School 
(NPS) data. 

The environmental Physics Group of NPS has conducted several 
research cruises over the past five years in both the Atlantic and 
East Pacific, at all times of the year, so that a considerable 


data base has been developed for a wide range of marine 


atmospheric conditions. Five of these cruises were conducted in 
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California coastal waters and the data needed to determine 
boundary layer mixing rates were collected. . 
These data have been gathered near enough to the coast to be 
representative of overwater cells in current air pollution models, 
and far enough at sea so that open ocean conditions prevail. The 
immediate areas of the three most populous centers, San Francisco, 
Los Angeles, and San Diego have been extensively investigated. A 
a wide range of meteorological conditions have been encountered, and 
operations were carried out around the clock. Thus, the available 
data base allows determinations of mixing rates and their 
dependence on location and conditions including diurnal 
: variations. In general, conditions over the water do not vary 
nearly as much as overland. In California coastal waters 
conditions are generally unstable with slightly stable encountered 


occasionally. 
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II. Instrumentation 
The studies of California coastal waters have all been carried 
out aboard the NPS ship R/V ACANIA, which is equipped to make 
multilevel mean and fluctuation measurements. Available sensor 
heights above mean sea level are 1, 4, 7, 14.5, 20.5 meters. For 
most of the data reported here, 4.2, 7.0, and 20.5 meters were 
used. Boundary layer properties were determined by three 
techniques: 
1. bulk aerodynamic (air-sea differences) 
2. inertial subrange turbulence 
3. mean profiles 
Rather than list all of the shipboard equipment here, we 
restrict the description to that which is pertinent to determining 
mixing rates. + The measurements needed and the equipment used 
were: 
Sea surface temperature (Ts) 
Air temperature (Tg) 
Humidity/Dew Point (Tp) 
Relative wind speed and direction (U) 
Temperature inversion height (24) 
Wind speed fluctuation (U') 
Air temperature fluctuation (T') 
Ty: 1) Hewlett Packard 2801A quartz thermometer ( 0.1°C) 
2) Barnes PRT-5 infrared thermometer ( 0.3°C) 


The HP sensor is floated and averages to about six ; 


inches below the surface. 
Tag: Same Hewlett Packard System (+0.2°C) 
Sensors installed in RM Young Gill aspirators. - 
The lower opening of the aspirator has been fitted 
with a radiation shield which improves its 
performance. 
Tp: General Eastern 1200 AP (Tp, +0.3°C) 
Dew point measured by cooled mirror technique 
Uo: MRI 1022 (+0.5 knt, +10°) 
The cups are low threshold so that 1 knt can be 
measured. Due to inaccuracies in ship speed true 
wind error is +1 knt, +15° at best. 
Z,: Aerovironment Model 200 acoustic sounder (+20m) 
Enclosure designed for shipboard allows good 
signal to noise when ship is in motion. 
U': TSI 1054B Hot Wire Anemometer 
Sensor is 1210 probe mounted with platinum film on 
quartz cylinder (60u) substrate. 
T': Sylvania 140 Thermosonde 
Sensor is TSI 1210 probe mounted with 2.54 
platinum wire. 
The equipment evolved with time and the equipment described above 
is the latest version used. 
Sensors were mounted on the R/V ACANIA so that the ship 
disturbs the sampled air as little as possible. The sensors at 1 


m height are placed on a bouy foreward of the ships bow (data 
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from this height is not used here). The 4 m and 7 m sensors were 
mounted forward of the ship on a mast placed directly on the bow. 
The 14.5 m level suffers the most ship influence and this data was 
only used under special circumstances. The 14.5 m and 20.5 m 
levels were located on a mast aproximately 15 ft. behind the bow. 
We attempt to obtain data only when the relative wind is within 
30° of the bow but this is not always possible. Even with the 
precautions of best sensor placement and good relative wind 
direction it has not been possible to obtain reliable wind 
profiles on the ship. 

Several methods of signal processing and data acquisition were 
used. This gives as much flexiblity as possible in choosing which 
of the three methods is used in the computations. Only the : 
fluctuation signals require significant processing. Two schemes 
were used: 1) spatial filtering and 2) frequency filtering. 
Spatial filtering was accomplished by placing two sensors 0.3 m 
apart and determining the difference in their responses. This 
method requires that matched sensors and processing electronics 
(both de and ac response) be used. Frequent checks on the senors 
were made to insure that environmental aging did not cause their 
responses to differ more than is tolerable. For frequency 
filtering, a single sensor was used and bandpass filtering, with 
lower and upper cutoff frequencies of 6 and 200 Hz, was imposed. 
The 0.3 m seperation and the 6 to 200 Hz bandpass both insure that 
only fluctuation components in the inertial subrange are utilized. 


It is necessary to restrict measurements to the inertial subrange 
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since ship motions introduce signal at lower frequencies, which 
would lead to incorrect results if these frequencies were used for 
direct flux estimates. 

After spatial or frequency filtering, the rms value of the 
fluctuation signal is obtained which is then used in subsequent 
calculations. It is very difficult to obtain matched sensors and 
to construct a difference bridge for the hot film sensors used 
for wind speed fluctuation measurements. Thus, only frequency 
filtering was used for wind speed fluctuations. Both types of 
filtering were used for temperature fluctuations, but the majority 
of the results are for the spatial filtering technique. 

The third method used to process fluctuation signals was 
spectral analysis. This can be done only for single sensors. The 
power spectral density of the inertial subrange signal was 
determined in-situ with a real time spectrum analyzer. The 
results were not one of the primary analysis tools, but were used 
to check the validity of the results from the other techniques. 
For example: 60 Hz pickup would be apparent by spectral analysis, 
but would increase the rms signal without the operator's 
Knowledge, and lead to erroneous results. Other equipment 
problems, and ship influence distortions can also be detected by 
Spectral analysis. Thus, spectra were produced on a frequent and 
regular basis throughout all cruises. Spectral analysis was 
performed by Nicolet 440B or Federal Scientific VA500 analyzers. 

Data acquisition was straightforward. Fluctuation signals 
were recorded on a Honeywell 5600 FM tape recorder. All mean and 


rms signals were acquired and recorded by a Hewlett Packard 3052A 
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data acquisition system controlled by an HP9825S computer. 
Complete cycling of the acquisition system through all signals 
took approximately 1.2 sec. Data was normally acquired for a 30 
minute period then averaged, so that each averaging period 
contains approximately 1200 samples. 

The computer performed in-situ calculations of meteorological 
parameters and recorded the data and results on magnetic tape. An 
immediate printout of all results was also produced. The 
immediate printout was very important for a successful operation 
as it allowed identification of portions of the system that were 


operating incorrectly. 


IIT. Details of Research Cruises 
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We report here results from five cruises: 


CEWCOM-76, 10/4-10/12, 1976; area covered was from Monterey Bay to 


San Diego with emphasis on the area south of Pt Conception. 
ARB, 7/19-7/26, 1977; area covered was Los Angeles air basin. 


CEWCOM-78, 5/14-5/23, 1978; coastal data taken but emphasis on at 


sea data, San Nicolas Island. 


MABLES-"C, 7/31~8/17, 1978; all data taken in San Francisco area 


from coast to 60 N mi at sea. 


Ctq, 6/2-6/8, 1979; Monterey Bay area, all data taken within 10 N 


mi of coast. 


Charts for these cruises follow so that the results can be 
correlated with the ships location with respect to the coast. 
Dates shown on the charts indicate the location at 0000, local 


time, except where the time is indicated. 


The chart shown for MABLES-WC covers only the time period from 
8/2 to 8/4. The remainder of the time was spent along latitude 
37°10'N at three stations located at longitudes: A=--122°40'W, 
B--123°15'W, C--123°50'W. The schedule for the three stations is 


ra shown in Table 1. 


Table 1. Ship schedule for MABLES-WC boundary layer study. 


Date/Time 0400 1000 1600 2200 
8/4 B Cc 
5 C B A A 
6 B C C B 
7 A A B C 
8 C B A A 
9 B Cc B B 
10 A B C C 
ame B A B Cc 
le C B A A 
13 B C C B 
14 A A B Cc 
1S C B A A 
16 B C C B 
17 A A 


Weather permitting, the ship sailed at full ahead between 


Cw tem te 


stations. It arrived at each station at least one hour before and 
left approximately one hour after the appointed time. The ship 

: was positioned downwind of the station at a distance such that it 
‘ would cross the station at the appointed time by sailing slow 
ahead into the wind. The charts for the CEWCOM-76 cruise do not 
mea. oe show the ships location from 10/9-1100 to 10/10-1200. ‘The only 

: . mixing rate data for which no position is shown occurred after 


0700 on 10/10. During this time the ship was immediately north 


eT 
ONT ee 


of Pt. Conception within 10 kmi of shore. 
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IV. Data Analysis 


The primary purpose of this work is to determine boundary 
layer mixing rates and to attempt to relate the rates in a 
consistent manner to local conditions. For each time period, the 
stability, Us, Ta, mixing rate, and mixing time were 
determined. We have found that the bulk and turbulence methods 
were the best to use for data analysis. Both methods were used to 
process these data but only the bulk results are reported. The 
turbulence results were used to check the validity of the output. 

In this section we describe the methods used to reduce the 
data. The spectral analysis technique will also be described for 
sake of completeness. The key to the analysis is the 
determination of hydrostatic stability and the resultant stability 
correction functions that must be used to determine final 
parameters. As will be seen, this is an iterative process that is 


best done by computer. 
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In what follows we use the notation: 


T 

Zz 

U 

q 

@ = T + .0098Z 
8) =@ + .61Tq 


temperature (K) 

height (m) 

mean wind speed (m/sec) 

water vapor mixing ratio 
potential temperature 

virtual potential temperature 
acceleration due to gravity 
von Karman's constant (0.35) 
wavenumber 

frequency 

turbulent diffusivity ratio (a, = 1) 
friction velocity (scaling velocity) 


mixing ratio scaling parameter 


( 


scaling temperature 

mixing rate 

roughness length 

Monin~Obukhov Tength 

inversion height 

atmospheric parameter (T, q or U) 
drag coefficient 

Z/L 

gradient stability function 


profile stability function 


stability function 


P(E) dissipation stability function 


€ turbulent kinetic dissipation rate 
Ct. structure function 
sy power spectral density 


We describe the stability with the Monin-Obukhov length 


2 
_ T Us 
Li Oye (1) 


The scaling parameters are related to the gradients by 


ae at P(e) (2) 


Current best evidence shows that transport of the scalars, heat and water 
vapor, obey the same relationship to their gradients. Thus, GY, * if and 


dp =a, = 1.35. The stability correction functions are:° 


q 
P(e) = (1 ~ 96)" B<0 © P(e) = (1 - 158) 
(3) 
= (1 + 6.42) —E>0 = (1 + 4.78) 


Bulk Method 


To obtain the stability the scaling wind speed and temperature are 
ilo k needed for use in Equation 1. They are obtained from bulk measurements 


using an integrated form of Equation 2. We integrate Equation 2 from the 


surface to a reference height Z, usually 10m. 


ae ee ee 


Z 


Z 
E 
‘ (8) 3 
0 


(4) 


_x z | 
oer [* 7 v4(e) 


where X, is the surface value. 


For analysis of data, it is most convenient to solve Equation 4 for X, 


and rewrite in terms of the drag coefficient. The neutral stability drag 


coefficient is given by 


5 aK 
Cw * wT zit, (5) 


and corrected for stability by 


ik byl 
Cy Qy [ - vy y/o a (6) 


Thus the scaling parameter is given by 


15 
Xe =C,(X, - X,) (7) 


As was stated above the key to the analysis is obtaining the stability. 


From Equation 1 we write & as 
g-f g Z Ox + 0.61 Tax (8) 
Us 


For the humidity correction term, we use T = 15°C and approximate 
0.61 T = 0.18 (q in gm/kg). Rewriting in drag coefficient form gives 
for neutral stability. 


4s 
: _«az Cw (T - T.) + 0.18 (q - q,) ' 
"oT OG u2 
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where we assume zero wind speed at the surface. Using Equations 6, 8, and 


9 gives 


x 2 
[1 - (EVE y/e] 


Ee, 


] ea Wy (EVEy/28 
The procedure used to analyze bulk data is as follows: 
1. Calculate q and q. from the measured Ty» T and relative humidity 
for Z = 10 m and assuming humidity - 100% at the surface. 


2. Use C,, = 1.3 x 1073 and Cun from the table below and the measure- 


TN 
ments to obtain Eo" 

3. Calculate by and by: 

4. Calculate & from Equation 10. 

5. Iterate steps 3 and 4 until the desired accuracy is obtained, 
giving &. 

6. L has been determined and U,, T,, and q, are obtained directly 
from Equation 7. 


4 
For Z = 10 m the wind drag coefficient is found from 


U_(m/sec) Cun Xx 10° 


0.3 - 2.2 1.08 y7!-5 


2.2 - 5.0 0.77 + 0.086U 
5 - 8 0.87 + 0.067U 
8 - 25 1.2 + 0.025U 


The profile stability functions are 


14x 
vrlé < 0) =2 un(—5") for x = (1 - 9&) 


Uy(E > 0) = -6.5¢ 
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vl <0) =2 an( 3X) + an(lexe —z—) - 2 tan “ly +5 
for x =(1-+ en 


v(E > 0) = - 4.78 


Turbulence Method 


(12) 


The turbulence method of data analysis uses the relationship between 


the structure function and the scaling parameter 


2_ 2 2-2/3 
Cc eX, 2" £6) 
with 
fee) = 4.9(1 - 76)"7/F gc Fle) = 4(1 + 0.5) 6 17/9) 
= 4.9(1 + 2.4¢2/3) E>0 =4(1 + 2.5¢2/3) 


It is more usual to use e rather than Cr They are related by 


The scaling velocity and the dissipation are related by 
3 
e€ = (Uy /eZ) P(E) 
with 


3/2 
PAE < 0) = (1 + 0.5/6/°/%) 


3/2 
PAE >O0)= (1+ 2.5¢2/3) 


(13) 


(14) 


(15) 


(16) 


(17) 
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Turbulence measurements in the inertial subrange yield - directly. 
The power spectral density, S.(k), of the inertial subrange portion of 


the turbulence is? 
S(k) = 0.25 es K79/3 (18) 


The spectral density is measured as a function of frequency. Using the 


frozen turbulence hypothesis k = 2nf/U and f S.(f) = k Sk) gives 


C2 = 4 s,(f)(2n/u)?/ ls (19) 


Measuring the squared mean difference signal from two sensors separated 


a distance d gives the structure function directly: 


C2 = Ixr) - x(n + ayyeya2/ (20) 


If the fluctuation signal from a single sensor is bandpass filtered at 
lower and upper wavenumbers kp and ky then the rms fluctuation signal 


Squared is 


u 4 
1 y2y 
(xX via! f S,(k) d k (21) 
kp 
Substituting Equation 18, integrating and using the frozen turbulence 
hypothesis gives 


-1 
2/3 sy, \o re “2/3 2/3 
c,? 2 ge) (Xing) (fp > fy od (22) 


Thus, a measurement of the rms signal yields cs directly. 


Mixing Rate 


The scaling parameters described above apply to, and specify the 
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state of the surface layer. The surface Tayer momentum, heat and water 
vapor fluxes can be determined from scaling parameters by 
2 


F_ = Us 


Fy = Us Ge (23) 


The Fi, in Equation 23 is the sensible heat flux. (The virtual heat flux 


and hence 8, is needed to calculate the mixing rate.) The depth over 
* 


which the surface scaling applies is approximately |L]. Above the surface 
layer (ignoring the transition is the well mixed layer region), surface 
scaling no longer applies and new scaling length, velocity, and temperature 


are needed. These are: 6 


Z, 


ws = C(a/T) Q, 2,113 
H, = Q,/» 
We assume that the boundary layer depth is defined by the inversion height, 


Z;. Q is the surface virtual heat flux. 


w, is the scaling velocity in the well mixed layer, and we assume 
that this velocity is directly related to the boundary layer mixing rate. 
Previous SF, tracer experiments performed by California Institute of 
Technology in Cooperation with NPS have shown that w, computed by the 
above method closely predicts the tracer experiment results. This is the 
currently available verification for using w, as the predicted mixing 


rate. 


The boundary layer depth, Z, used in Equs. 24 cannot always 
be taken as the inversion height. During times when there is a 
stratus layer, the lower edge of the stratus forms a boundary 
above which condensation must be taken into account in the heat 
balance. For such conditions Equs. 24 apply to the volume of air 
below the stratus, and this height should be used for Z,. If 
the inversion height is used, the calculated mixing rate will be 
somewhat in error. This will not be serious for most conditions 


since the height appears to the one third power in the expression. 
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RESULTS 

There is no accepted, nor simple, manner to present the 
results of a study of this type. The parameter of interest, the 
mixing rate, depends on the full range of atmospheric parameters, 
wind, heat flux, stability, boundary layer height, etc. These 
parameters are influenced by proximity to land, which can also 
lead to diurnal variations due to the solar cycle. Since all 
parameters can change with time, and their variations are not 
simply related, one cannot relate the mixing rate to one or two 
basic parameters. Of course given a complete set of 
meteorological measurements the mixing rate can be calculated from 
Equation 24. 

Our purpose here is twofold: 1) provide a simplified method 
to calculate the mixing rate from locally measured meteorological 
parameters, 2) catalog average values of the mixing rate that can 
be used to improve current air pollution models. The average 
values presented will be appropriate to location, time, and 
generalized local conditions. 

Tabular data and results are presented in appendices A and 
B. Appendix A is the basic meteorological data: wind speed, 
relative humidity, air and sea temperature, and inversion height. 
The air sea temperature difference is included since it and the 
wind speed are the determining factors for the stability and heat 
flux. The final entry is the heat flux, which together with the 


inversion height, gives the mixing rate (Equation 24). Appendix 


B consists mainly of calculated parameters. The wind speed and 
direction, and the inversion height are included as references. 
The wind direction is especially useful here since it can be used 
to differentiate local (land and sea breeze) circulations. The 
tables also include the stability (Z/L), the scaling parameters 
Ux, Tx, the mixing rate wy, and the total mixing time. The 

total mixing time is the time it takes to move a parcel of air 
from the surface to the top of the marine layer (or vice-versa). 

Note that there are many entries missing from the tables. 
This is due to one sensor failing during a time period or to 
stable conditions. The theory is invalid for stable conditions, 
which occur a very small percentage of the time. The valid data 
for the time period is included in order to maintain as complete a 
record as possible. 

Figures 5~9 show all of the mixing rate results, where the 
rate is plotted as a function of time for the full period of each 
cruise. All plots except MABLES-WC show these results for 
one-half hour averaging periods. The shaded areas at the bottom 
of the graphs show times when the ship was near shore. Since the 
Ctq and ARB cruises were mostly near shore, no shading is 
shown. The MABLES-WC data appears in a special format which will 
be discussed seperately. 7 

Several general results are immediately apparent from a 
cursory examination of these plotted results: 

1) Value range from 0.2 to 0.8 m/sec 

2) Low values generally occur near land 


3) Fluctuations from one one-half hour period to the next as 
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large as 0.2 m/sec are to be expected. 

The fluctuations are attributed to both changes in wind speed 
and air-sea temperature difference. In many cases the ship was 
underway while the data was being taken, which accounts variations 
in observed sea surface temperature, and hence air-sea temperature 
difference, over short periods of time. When the ship was at a 
stationary location, large changes in mixing rate due to wind 
speed fluctuations were observed. This was especially evident 
near the shore. 

In general, conditions along the California coast were found 


to be quite consistent both spatially and temporally. We now 


‘summarize the principle features of the results from each cruise 


in order to illustrate the important parameters. 


CEWCOM-76 (Figure 5): The ship operated in the area extending 
from about 100 N mi north of Pt. Conception to San Diego for the 
period covered in this report. Near shore values of ws were 
approximately 0.3 m/sec except for 10/6 when the value was ‘0.5 
m/sec due to a fairly large air-sea temperature difference 

(AT ~ 2°C). The low values on 10/5 occured when wind speed and 
air-sea temperature difference were simultaneously low. The rapid 
decrease of wa in the afternoon of 10/7 was associated with the 
lowering of the inversion base as the ship neared shore (refer to 
Equation 24). The larger values on 10/9 were associated with a 


large AT v3°C. 


ARB (Figure 6): The ship operated in the Los Angeles to San Diego 
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area, spending most of the time near shore. The largest ws of 
approximately 0.6 m/sec was observed when the ship was away from 
the shore on 7/19 and was associated with large AT values. Near 
shore values averaged 70.3 m/sec. The very low values (< 0.2 
m/sec) occured when AT was small, except on 7/22 when low winds 
resulted in low heat flux anduwsx. The very rapid decrease on 

the morning of the 19th was also associated with a decreasing wind 


speed in Santa Monica Bay. 


CEWCOM-78 (Figure 7): This ship operated in the vicinity of San 
Nicolas Island, the channel islands and the Santa Barbara 
shoreline. Again, the lowest values for ws were obtained when 

the ship was near land. Low values, ~0.15 m/sec, on 5/14 occured 
in near neutral conditions. Extremely high winds were encountered 
on 5/15 (up to 60 knts, when the anehon shad broke) and ws 
increased rapidly. Large values, v0.7 m/sec, from 5/19 to 5/22 
were associated with large AT values. AT, and hence ws, 


decreased as the ship moved to the mainland on 5/22. 


Ctg (Figure 7): All data was taken near shore in Monterey 
Bay. The largest wx 0.75 m/sec occured for a large AT v2.5°C 
and moderate winds. Lowest values and large fluctuations occured 


on 10/7 when conditions were near neutral. 


MABLES~WC (Figure 8): The data shown is for the period when the 
Ship was operating at stations A, B and C as described in the 


previous section. The data is averaged for the full time the ship 
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was at either station A or C, and both sets of results are plotted 
as functions of time. This allows an immediate comparison to be 
made between offshore (C) and near shore (A) locations. It is 
immediately apparent that offshore values of ws are generally 
higher than near shore values. This is due to generally higher 
offshore wind speeds. The only time when ws was higher near 
shore (on 8/14) was when there was a very high AT ~~ 3°C. Along 
the California coast offshore winds are generally higher than 


winds near the shore. 


¢ 


We now summarize these results in a form appropriate for easy 
use in air pollution models. Table II presents average values of 
Wy appropriate for various conditions and locations. These 
values can be used in place of a value calculated from measured 
conditions. Accuracy can not be expected to be any greater than 
50% and an average avlue could be a factor of two different than 
the actual value. | 

TABLE II 


Average mixing rates for various locations and conditions 


Conditions Rate (m/sec) 
Open Ocean 0.7 
Within 10 N mi of coast 

24 hour average 0.4 
Night 0.5 
Afternoon 0.6 
Land-Sea breeze changeover 0.3 
Strong northwesterly, any location 0.7 


Low wind and air-sea temperature difference 0.2 


- MABLES-WC > 
(AUG-1978) 


It must be reemphasized that the above recommended values for 
the mixing rate must be used with care. There are several 
commonly encountered conditions that will cause the rate to be 
different than that given in the table. In particular the table 
does not take the air-sea temperature difference, nor the 
magnitude of the wind, into account. We have observed wa to be 
as large as 0.9 m/sec with a large temperature difference in 
moderate winds. During the winter there is often little diurnal 
variation in the wind and a value wy % 0.5 m/sec for all times 
is reasonable. 

Of course, the most accurate means of determining the mixing 
rate is to use Equation 24 with measured mean meteorological data. 
Full utilization of this method requires the use of the correct 
wind dependent drag coefficient and an iterative scheme to 
determine the stability. A simplified calculation scheme can be 
used that is fairly accurate. We now outline this method and the 
final equation that replaces Equation 24. 

Approximate drag coefficients of €, = C, = 1.3 x 10-3 
can be used with the bulk method. These values are substituted 
into Equation 7 to give approximate scaling parameters Ux, Ts 
ax. The virtual heat flux can then be obtained directly from 


Equation 23 and used to calculate the mixing rate. The result is 
/3 
Wy, = 0.035 [2,19 [(T9 - Ts) + 0-18(aqyq - Ig) ] (25) 


where Uj0, Tio and qjo refer to the values measured 


at 10 meters above sea level. 


SO Ta neo RT TORTIE TEER ETO Cu OeETS 


Measurements may be made at a height other than 10 meters. 
The drag coefficients can be easily adjusted to any reference 


height. Using Equation 5 we have 
ef = ef (in 10/26)/ (1n 2/20) (26) 


where Zo, the roughness length, is approximated by 6 x 1074 

m for wind and 2 x 1072 m for temperature and water vapor. 

Note that corrections to Equation 25 must be applied for both the 
wind, and for temperature and humidity. The resulting correction 


forws, is 


w,(Z) = w,(Z = 10) 5.03[(1nz + 7.42) (Anz + 10.82) ]7 “2 (27) 


If a determination of sea surface temperature and inversion 
height is available, then the mean air parameters can be 
determined on shore in a region where acceleration effects are 
small. These measured parameters can be used in Equation 25 to 
give the mixing rate with an expected error of no greater than 
25%. The expected error is due to ignoring the dependence of the 
wind drag coefficient on wind speed and the slight dependence of 
all drag coefficients on hydrostatic stability in the unstable 


regime. 
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Appendix A 


Meteorological data: The data are arranged in chronological order 
for each of the five cruises. Included are wind speed, relative 
humidity, air temperature, sea surface temperature, air-sea 


temperature difference, and the calculated virtual heat flux. 
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Appendix 


Calenlated results: Che results are 


order ror each of the five crufses. 


direction, stabllity (2/L), the sealing parameters Us and 7 
inverston heteht, mixing rate, and mixing 


appears after the mixing time, it means that the relative wind was 


more than 30° off the bow. 
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Included are 
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